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The particle spectra and Hanbury-Brown Twiss (HBT) radius ol Au+Au collisions at RHIC 
energy are investigated by a hydrodynamical expanding source with both shear and bulk viscosities 
(£). With a large width of the ratio of £ to entropy density s, both the particle transverse momentum 
spectra and the ratio R out / R S ide of HBT radii in the direction of the total transverse momentum 
of detected two particles (Rout) and perpendicular to both this direction and the beam direction 
(Rside) become a little steeper. 
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I. INTRODUCTION 

Ideal hydrodynamics has successfully explained the ex- 
perimental data of elliptic flow coefficient V2(pr) as well 
as the transverse momentum spectra of hadrons in central 
and semi-peripheral collisions up to px = 1.5-2 GeV/c 
and also predicted the hard photon production. I 1 ' 2 1 How- 
ever the almost perfect ideal fluid dynamical descrip- 
tion of the experimental data gradually breaks down 
in the regions of more peripheral collisions and/or high 
transverse momenta (pr), and it also fails to repro- 
duce the HBT radius J 3,4 ' Actually ideal fluid is only an 
approximation, I 5 1 people need to investigate the viscous 
hydrodynamics. The so-called the first order relativistic 
viscous hydrodynamics was firstly formulated by Eckart 
and a variation by Landau and Lifshitz, and many cal- 
culations have been done so far J 6-8 ' But there were two 
problems in their approaches: one is the dissipative fluc- 
tuations which may propagate at a speed larger than the 
speed of light and thus lead to the causality problem; 
the other is the solutions which may develop instabili- 
ties. The second order formalism which was developed 
about thirty years ago by Israel and Stewart avoided 
these difficulties. I 9 ! But due to the complication of the 
evolution equations, numerical codes have not been de- 
veloped to solve the Israel-Stewart equations until recent 
years. I 10-16 ! When compared to the experimental data, 
people only considered the shear viscosity and found that 
the particle spectra and elliptic flow became more con- 
sistent in larger px range. But the HBT radius are still 
not consistent with the experimental data. I 12 ! 

Although one expects that the bulk viscosity is to be 
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negligible for temperatures far away from the phase tran- 
sition temperature T c , large effects are expected near 



T c } 17 ~ 21 ' As suggested in |22|, the HBT results may come 
close to the experimental data if bulk viscosity is taken 
into account. Based on this argument, we want to check 
if this argument is valid in the present work. To this end, 
we use the algorithm which has been used in [12j to con- 
sider the effect of the bulk viscosity for central Au+Au 
collision at RHIC energies. 



II. CAUSAL VISCOUS HYDRODYNAMICS 
AND RESULTS 



The general hydrodynamic equations arise from the 
local conservation of energy and momentum! 1 ] 



d^"{x) = 0, 



(1) 



where the energy-momentum tensor without heat con- 
duction is decomposed in the form! 14 ! 



(2) 



Here e and p are the local energy density and thermal 
equilibrium pressure, and is the 4-velocity of the en- 
ergy flow which obeys u^u^ — 1. n is the bulk viscous 
pressure, A Miy = g^ v — u^u v is transverse to the flow ve- 
locity, that is A^iiv — 0. ir^ is the traceless shear vis- 
cous pressure tensor. With Eq. ((T|) and Eq. ([2]), we can 
get the evolution equations of the energy density and the 
4-velocity of the energy flow 

{e+p + U)Du' i = V*(p + n) 

- AJJV^tt'" 7 + ir^Du v , (3) 



De = -( e + p + n)V M u" 



(4) 



where D is the time derivative in the local fluid rest frame 
and fulfills D = u^d^. The angular bracket notation is 
defined by (Wu u ) = ^ ( V^+Vu") - ±(V-u)A^. In 
the Israel-Stewart approach the kinetic evolution equa- 
tions of the bulk pressure n and the traceless shear vis- 
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cous tensor are I 14 ' 

DU = — (n + cv-u), (5) 

Dtt^ = - — {^"-2rf7^u v) ) 

-(u^vr^ + u l, n^ a )Du a , (6) 

where rj and £ denotes the bulk and shear viscous coeffi- 
cients, respectively, and rn are the relaxation times for 
the bulk pressure and the shear viscous pressure tensor, 

respectively. They can be related to 77 and ( as follows 

[10] 

7v=2ry/3 2 , m = CA>- (7) 

For a massless Boltzmann gas ~ so we can get 
7V sy which translates to w ^^l 10 ' 11 ! where s is 
entropy density and this is the value we use in this paper. 
In this case the relaxation and coupling coefficients for 
bulk pressure vanish, /3 — > 00 and thus there are no ef- 
fects of the bulk viscosity. ^ But as we already know the 
bulk viscous coefficient £ becomes large near the quark- 
hadron phase transition temperature T C J 23 1 So the be- 
havior of the bulk viscous pressure must change and /3q 
does not tend to be infinite any more. However, as to our 
knowledge, there are no quantitative results for /3o that 
are relevant for applications to nuclear collisions (an at- 
tempt has been made in [To|). In this paper we will try 
two different /3q in our calculations which correspond to 
two different Tn , namely, rn = ~ and rn = § |ft , respec- 
tively. Of course people could try smaller /3 . But with a 
smaller /3 , there will be casuality problem and this may 
arise instability. I 24 ' 

As in ideal hydrodynamics we need initial conditions 
at To, equation of state of the bulk matter and freeze- 
out condition to obtain the results of these evolution 
equations. Generally there are two models for the ini- 
tialization of the hydrodynamics: Glauber-type model 
and Color Glass Condensate-type model. Here we use 
the Glauber model and the energy density at the ini- 
tial time To is parameterized by the number density of 
wounded nucleons e ~ hwn,' 12,25 ' where n\vN = 27a [1 — 
(1 - ^) A ] and T A = dWU + exp[(V^T^ - 
Ro)/x]}- For gold nuclei, A = 197, R = 6.4 fm, X = 0.54 
fm, the nucleon-nucleon cross section a at ^Jsnn = 200 
GeV is assumed to be 40 mb and po is chosen to fulfill 
2ir J rdrTA = A. The EOS we use here is calculated by 
Laine and Schroder which corresponds to a phase tran- 
sition temperature T c m 190 ± 20 MeV. [261 In their cal- 
culation they considered a hadron resonance gas at low 
temperatures and used high order weak-coupling QCD 
result at high temperatures. And the freeze-out mech- 
anism is the Cooper-Frye prescription that is particles 
freeze out at a single temperature T/J 27 l 
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Fig.l. Two linearly changed with different width and 
height. See texts for details. 

Unlike ideal hydrodynamics, in causal viscous hydro- 
dynamics we need to give the initial condition for t:^ u and 
n as well as the values of rj/s and C/sJ 10-12 ' 14 ' Here we 
choose tt^ = and n = at the initial time which min- 
imize the effects of viscosity. There have been many cal- 
culations for r\j s in which there is a conjectured minimal 
bound for r]/s in AdS/CFT that 77/s > 1/4tt « 0.08, 1 28 l 
and we can see in [29| that r\j s does not change a lot for 
different temperatures. In this paper we want to see the 
influence of the bulk viscosity, so we use the lower bound 
of 77/s during the whole evolution of the fluid. In [2l[ 
the value of £/s reaches the maximum at T c and changes 
linearly to a very small value when the temperatu re g oes 
away from T c and the width is small. However, in [T9II20II 
the width is large. In this paper, we will take two linearly 
changing £/s with different width and height as shown 
in Fig.l in our calculations. 

For central Au+Au collisions at RHIC, the geometries 
which are longitudinally expanding are space-time rapid- 
ity independent and have radial symmetry. For these ge- 
ometries it is convenient to work in the co-moving and ra- 
dial coordinates t, r, <j), rj with the relations t = y/t 2 — z 2 , 
r 2 = x 2 + y 2 , tan0 = y/x and 77 = atanh(z//j). 

The only non-vanishing fluid velocity components are 
then u T and u r with the relation u T = v/l + (u r ) 2 . To 
solve Israel- Stewart hydrodynamic equations which have 
this kind of geometry we use the one dimensional algo- 
rithm used in (T^ |. 

As the procedure that adopted in the ideal 
hydrodynamics, I 30 ' the initial central temperature To and 
the freeze-out temperature Tf are chosen to make both 
the normalization and the slope of the resulting pion 
spectrum in reasonable agreement with the experimen- 
tal data and we get T Q = 0.34 MeV and Tf = 0.17 MeV. 
The freeze-out temperature is larger than that typically 
used in ideal hydrodynamics. In ideal hydrodynamics in- 
teractions of the system can keep the system in thermal 
equilibrium, but the interactions are not so efficient in 
viscous hydrodynamics which allows for departures from 
equilibrium. So an earlier freeze-out and thus a higher Tf 
is to be expected in viscous hydrodynamics. In [131 ] there 
are other parameters that can fit both the pion spectrum 
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and the HBT radii, but either the r\j s is too large or the 
freeze-out temperature is unreasonable large. 

Fig. 2 shows the velocity and temperature profiles for 
different times of the evolution of the fluid using the 
parameters after fitting the pion spectra. The left and 
right panels are for relaxation time values rn = % an d 
Tn = - ^r, respectively. In each panel, we calculate three 
cases: (a) rj/s = 0.08 and without bulk viscosity, (b) 
rj/s = 0.08, (^) m ax — 0.3 and the later changes lin- 
early to (§) mm = 0.005 at 0.98T C and 1.02T C respec- 
tively, corresponding to the solid line in Fig.l, and (c) 
■q/s = 0.08, (^) m ax — 0.1 and the later changes linearly 
to (^)min = 0.005 at 0.9T C and 1.2T C respectively, cor- 
responding to the dashed line in Fig.l. We can see that 
if bulk viscosity is taken into account, at earlier time 
the velocity changes near T c , and in case (c) the veloc- 
ity changes more rapidly than case (b). This is due to 
the narrower width and larger height of £/s. And with a 
much larger peak there may be non-periodic oscillation 
in velocity and temperature. So the method of hydrody- 
namics may not be suitable and the source may clusterize 
to small fragments.! 31 ] At later time, the effect of bulk 
viscosity extends to other place where the temperature 
is far away T c . We can see that the evolution of the fluid 
is sensitive to the width of £/s rather than its height. In 
both case (b) and (c), the fluid tends to be retarded and 
the effect is stronger with larger width. We can also see 
that the effect of bulk viscosity is a little stronger with a 
smaller relaxation time. 
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Fig. 2. Temperature and velocity profiles for different conditions 
(a), (b) and (c) with rn = ^ (left panels) and tjt = 
(right panels) of the relaxation time values at different times 
T n = 2nfm/c. See texts for details. 

We plot the pion, kaon and proton spectra together 
with the RHIC STAR and PHENIX experimental data 
for the most central 5% of Au 4- Au collisions at ^/snn = 



200 GeV in Fig.3. [32 ' 331 It shows that case (a) and (b) are 
almost the same and case (c) become steeper which indi- 
cates smaller transverse flow. From this we can conclude 
that the effect of the width is more important than the 
height of C/s in hydrodynamics. 
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Fig. 3. Comparison of the calculated spectra of it, K (scaled by 
0.1) and P (scaled by 0.01) for different conditions (a), (b) and 
(c) with the STAR and PHENIX data. The upper and bottom 



panel corresponds to Tn = 
time values, respectively. 



and rn 
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Finally we also show the viscous hydrodynamic results 
for the HBT radii together with experimental data. The 
results are shown in Fig. 4 where the ratios R ou t I R side are 
for different relaxation time values and different cases. I 34 l 
The radius R ou t represents the radius in the direction of 
the total transverse momentum of detected two parti- 
cles and Rside corresponds to the radius perpendicular 
to both the total transverse momentum direction and 
the beam direction. There are some quantitative effects 
for HBT radius' ratio in different pr range. Essentially 
Rout/ Rside becomes steeper versus the transverse mo- 
mentum when the bulk viscosity is taken into account, 
which is consistent with the steeper transverse momen- 
tum spectra observed in Fig. 3. Nevertheless we still did 
not see the bulk viscosity can ch ang e the ratio substan- 
tially. This also can be seen in [351 ] . The reason could 
be that the thermodynamical quantities change only a 
little at Tf for different cases, and the freeze-out surface 
is almost the same. 



III. SUMMARY AND DISCUSSIONS 

In summary, we have solved the causal viscous hy- 
drodynamic equations for central Au+Au collisions at 
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^/snn = 200 GeV using a simple numerical code. Using 
the parameters that fitted the pion spectra, we 
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Fig. 4. Calculated ratio Rout / R a id<t °f HBT radii for different 
conditions (a), (b) and (c) together with STAR data. The upper 
and bottom panel corresponds to Tn = -I and rn = - of 
the relaxation time values, respectively. 



calculated the temperature and velocity profiles, the 
pion, kaon and proton spectra as well as the HBT results 
for different relaxation time values and different bulk vis- 
cosity cases. We find that the fluid evolves slowly when 
the bulk viscosity is taken into account. The particle 
spectra become steeper which indicates smaller trans- 
verse flow when the width of £/s becomes larger. We 
can also see that the HBT results are still far away from 
the experimental data when the bulk viscosity is taken 
into account in our approach. 

There are other problems that have to be considered 
to obtain a proper description of the experimental HBT 
radii. Firstly, we should investigate the influence of 
different EOS. In [22J the author use different EOS to 
calculate the causal viscous hydrodynamics and find the 
behavior of R ut/Rside for different £/s changes sub- 
stantially. Secondly, larger £/s may develop instabilities 
in the mixed phase, and the source may be clusterized to 
droplets and this freeze-out mechanism is different from 
Cooper- Frye formalism J 31 ' 36 ' A granular source method 
was claimed that it can explain the HBT puzzle J 37 ' 
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